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Abstract 

A conjecture due to Bohigas, Giannoni and Schmit (BGS), stating 
that the energy level correlations of quantum chaotic systems gener- 
ically obey the laws of random matrix theory, is given a precise for- 
mulation for quantized symplectic maps. No statement is made about 
any individual quantum map. Rather, a few-parameter ensemble of 
maps is considered, such that the deterministic map is composed with 
a diffusion operator on average. The ensemble is a "quantum" one, 
which is to say that the diffusion operator contracts to the identity 
in the classical limit. It is argued that the BGS conjecture is true 
on average over such an ensemble, provided that the classical map is 
mixing. The method used is closely related to the super symmetric 
formalism of Andreev et al for chaotic Hamiltonian systems. 



1 Introduction 

One of the goals of semiclassical analysis is to characterize the energy level 
correlations of quantum mechanical systems in the limit h ^ 0. To gain an 
understanding of these correlations, physicists have traditionally used semi- 
classical trace formulas, relating the density of energy levels to the periodic 
orbits of the classical dynamics (see e.g. the contributions by Bogomolny, 
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Keating, and Smilansky to this volume). In brief, the quahtative picture 
that has emerged is that integrable systems have Poisson statistics, whereas 
for a generic chaotic system one expects random matrix (or Wigner-Dyson) 
statistics. This picture is supported by a large body of numerical evidence, 
and applies equally to Hamiltonian systems and symplectic maps. Bohigas, 
Giannoni, and Schmit (BGS) are credited for having put forth the random 
matrix conjecture for quantum chaotic systems. The principal development 
toward an analytical theory is due to Berry who augmented Gutzwiller's 
trace formula with sum rule arguments to determine the form factor for times 
much shorter and much longer than the Heisenberg time. Berry's work has 
recently been refined by Bogomolny and Keating [Q]. 

In these lectures, I will describe a supersymmetric formalism that offers 
an alternative to the trace formula approach and promises rigorous results for 
chaotic symplectic maps. The formalism bears much similarity to that of An- 
dreev et al ||^ for chaotic Hamiltonian systems, reviewed at this workshop by 
Ben Simons. Only the simplest case of maps without unitary or anti-unitary 
symmetries will be treated here. In particular, time reversal invariance will 
be assumed to be broken. The starting point of our theory is an expression 
for the two-level correlation function which involves two ratios of quantum 
spectral determinants (Sect. The determinants are written as traces over 
a Fock space of fermions, and the inverse determinants as traces over a Fock 
space of bosons (Sect. H). By the introduction of generalized coherent states, 
the two-level correlation function can then be expressed as a Berezin inte- 
gral over the Riemannian symmetric superspace of type y4III|y4III (Sect. 
A quick review of the basic mathematics underlying the Berezin integral is 
given in Sect. |[ 

The BGS conjecture is not a precise statement, as the meaning of the 
word "generic" is left undefined. (The qualification "generic" serves to ex- 
clude some prominent counterexamples, such as arithmetic billiards and the 
cat maps, which are paradigms of classical chaos but fail to obey random 
matrix statistics when quantized in the canonical way.) It turns out that in 
our formalism it is impossible to prove the conjecture in the sense originally 
intended, i.e. as a statement about an individual quantum system. The tech- 
nical obstacle that prevents us from making progress for individual systems 
is the nonexistence of the semiclassical limit (Sect. ^ for the Berezin integral 
representation of the two-level correlation function, or any other correlation 
function. To enforce such a limit one needs to impose some sort of coarse 
graining, or ultraviolet regulator, on the Berezin integral. The way to im- 
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plement regularization is to compose the map with a stochastic Hamiltonian 
flow and average (Sect. ||). 

From experience with the supersymmetry formahsm apphed to disordered 
electron systems 0] it is known that the so-called zero mode approximation 
to the Berezin integral gives the random matrix answer. The strategy there- 
fore is to establish sufficient conditions for this approximation, which is of the 
saddle point type, to be justified. The key issue now becomes the stability 
of the manifold of saddle points. Without ensemble averaging, the Hessian 
of the saddle point manifold is not positive, and the saddle point approxi- 
mation is false. Averaging over the stochastic Hamiltonian fiow with "time" 
parameter e stabilizes the saddle point manifold and leads to a good asymp- 
totic expansion (Sect. |^). If the classical map is mixing we argue from power 
counting that stability can achieved with a time parameter e that vanishes 
with h as where < a < 1. Thus we propose that the BGS conjecture 
holds on average over an ensemble of quantum chaotic maps all of which have 
the same classical limit. 

My original intention was to finish with a discussion of the finite h correc- 
tions to the universal random matrix result. These are a very timely issue, as 
the ballistic nonlinear a model calculation of Andreev et al disagrees with the 
result of Bogomolny and Keating derived from the diagonal approximation. 
(The difference lies in the weights given to the contributions from periodic 
orbits with repetition number r > 2.) There exist strong indications that 
the result of Bogomolny and Keating is the correct one. I believe that the 
resolution of the discrepancy will teach us something profound about the bal- 
listic nonlinear a model, but this definitely exceeds the scope of the present 
lectures and must be left for a future publication. 

2 Preliminaries 

We begin with some basic definitions. 

Let M, called the phase space, be a 2(i-dimensional compact manifold 
with a symplectic structure u, i.e. an antisymmetric second rank tensor 
that is expressed in local coordinates g^, pk hy to = Y.k=i^Pk /\ dg^. Via 
u, every function / : M ^ M is associated with a Hamiltonian vector field 

= J2k=i ~ Mk'Sk)' '^^^ volume element, or Liouville measure, on 

M is denoted hj dx := u^'^. A map x : M — >■ M is called symplectic (or 
area preserving for d = 1) if it preserves u and, consequently, the Liouville 
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measure. The map x ^-cts on functions / : M — M by (x*/)(a;) = 

The dual of its inverse, x~^* i is called the Frobenius-Perron operator. A 

little computation, 

= / dy5{x'^{x)-y)f{y)= dy S{x - xiy))f{y) , 

JM JM 

shows that the Frobenius-Perron operator has the integral kernel 

{x'^*){x,y) = S{x - x{y)) , 

which identifies x^^* propagator (for one time step) of the classical 

dynamics obtained by iterating x- 

Quantization turns every smooth function / : M — R into a self-adjoint 
operator QnW) = -^a^ ^ Hilbert space Hn of finite dimension = 
{2TTh)~'^ Jj^^ dx. The quantization of a smooth phase space function is called 
a pseudodifferential operator. The quantum analog of the classical map 
X is a unitary operator Un '■ Ti-N — >■ Htv, whose eigenvalues we denote 
by e*^" {n = 1,...,N). The 6nS are sometimes called quasienergies. Let 
'D]\r = Hn ® '^jv — End(7iAr) be the space of linear operators on Hn- The 
unitary operator U n acts on Vn by conjugation. More precisely, the quantum 
time evolution (by one time step) of an element Fn € Vn is given by 

Ad{UN)-'FN = U^'FnUn . 

(The notation Ad{U) for conjugation by U is taken from Lie algebra theory.) 
In the classical limit h —>■ 0, Ad(f/7v)~^ acting on a pseudodifferential operator 
Fn = Qn^I) approaches the Frobenius-Perron operator: 

M{Un)-'Fn ^ Qiv(x-'V) • 

Examples illustrating these general definitions are furnished by the cat 
maps 0. The phase space in this case is the two-torus T^, coordinatized by 
a pair of canonically conjugate position and momentum functions q and p, 
which are defined modulo 1. For a concrete example, consider the cat map 
X : T^ — >• T^ that acts on q and p as 

1° X = q + '2p (mod 1) , 
pox = 2q + 5p (mod 1) , 
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which entails a Probenius-Perron operator expressed by 



X ^* = exp (—2qd/dp) exp (—2pd/dq) 

on differentiable functions. [Note that exp {—d/dq) = exp (—d/dp) = 1, so 
that exp{—2pd/dq) and cxp{—2qd/dp) arc globally well-defined, although 
pd/dq and qd/dp are not.] The locally defined vector fields = 2qd/dp 
and 'Eg = 2pd/dq derive from locally defined functions f — —(f and g = p^, 
and quantization of x yields a unitary operator 

Un = exp (iP^/hj exp (^-iQlf/hj , 

where — 2t^N Jq dp Jq dq — 2t^N. In a basis where exp (27riPjv) is diag- 
onal with the eigenvalues being the N^^ roots of unity exp {2Trin/N) {n — 
1,2, . . . , N), Un has the following matrix 

{UN)nn' = -7=exp — [n +{n-n) J . 



3 BGS conjecture formulated precisely 

There exist many statistical measures of level correlations. The one we will 
focus on here is the so-called pair correlation function, defined on a test 
function / : U(l) ^ C by 

rra,n=l q 

If we take the liberty of choosing for / the Dirac ^-distribution centered at 
e*^, 5^(e*^) = {2ti)~^ J2iez^^''^^~'^\ S^^ the two-level correlation function 



cy OO 

R2{9) = C{de) = ^ Ecos (W) TrU' 

1=1 



For our purposes, it is technically convenient to consider in addition to R2 
the following correlator of determinants: 



^(7+0, 7+1; 7-0, 7-1) = J ^ 





d9 Det(l - 7+ie*''f/jv)Dct(l - 7_ie-'^f/jv ; 
2^ Det(l - 7+oe^^t/iv)Det(l - 7_oe-*^t/^^) ' 



where 7±o,7±i are complex parameters. Note that the two-level correlation 
function can be extracted from Q by taking two derivatives: 



In Sects. I will sketch the derivation of a formula for fl that underlies 
the super symmetric method for quantized maps, and I will go over some of 
the basic mathematics needed. Before doing so, it is proper to motivate the 
reader by giving an outline of the results that can be obtained. 

The first and very important message is that our method gives no answer 
at all for the pair correlation function, or any other correlation function, of 
an individual quantized map. (This is to be contrasted with what has been 
claimed by Andreev et al.) For reasons that will be explained later on, we 
can make a statement only about the expected value of the correlations for 
an ensemble of quantized maps over a well-chosen probability space. 

Let me therefore digress and describe this ensemble. We pick a finite 
number s of Hamiltonian functions Xk '■ M ^ R, with associated Hamiltonian 
vector fields H^, and Hamiltonian operators X^^n '■ 'Hn T^n- The latter 
we multiply by uncorrelated Gaussian random variables C,k with zero mean 
and ?i-dependent variances (^|) = 2e{h), and we consider instead of a single 
Un an s-parameter family of unitary operators 



The choice of Hamiltonians Xk is constrained by the requirement that the 
weighted sum of squares J2k=i '^1 be an elliptic operator. 

The good news is that we can get an answer for the expected correlation 
functions (Q) and (-R2) by taking e{h) = const x ff where a is a number 
between and 1 (and we may choose a = 1/2, say). This dependence of 
the variances on h means that, roughly speaking, all members of the ensem- 
ble have the same classical limit. In other words, the ensemble we use is 
"quantum" , and we could say that what we are averaging over are different 
quantizations of the same classical map. 

Now, let the classical map be mixing, and arrange for the quantized map 
to have no unitary or anti-unitary symmetries. We will argue that the ex- 
pected correlation functions in this case have the following semiclassical lim- 
its: 
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(c+o - c_o)(c+i - C_i) ^ ^ 

and 

lim {R2{2nx/N)) = - sin^ {nx)/{nxf . 

N—*oo 

To put this result into context, we observe that our expressions for (Q) and 
{R2) coincide with those for a random unitary matrix in the large-A?^ hmit. 
What wc propose, then, is a precise formulation of the BGS conjecture: 
random matrix theory applies to chaotic maps on average over a quantum 
ensemble of the kind specified. More specifically, although there exist coun- 
terexamples such as the cat maps, these are a set of measure zero in parameter 
space, and the random matrix result does obtain on averaging. 



4 Supersymmetric formula for maps (I) 

In what follows we shall construct a novel and useful integral representa- 
tion of the correlator fl. Our first step will be to express the determinants 
Det(l — 7t/Ar)^^ as traces over Fock space. To that end, let J^n be the 2''^- 
dimensional fermionic Fock space that is generated by acting with creation 
operators fl {n = 1,2, . . . , N) on a particle vacuum |0). The operators 
fl and their adjoints satisfy the canonical anticommutation relation for 
fermions: + /„,/t = s„^,. 

Now we claim that determinants such as those appearing in the numerator 
of the defining expression for Q can be written as Fock space traces: 

N 

Det(l-C/^)=IV^^(-l)^-exp ^ /t(ln C/^)„„,/„, 

n,n'=l 

where Np — flf^ is the fermion number. To prove that claim, we may 
assume the matrix Un to be diagonal as both sides of the equation are in- 
variant under a change of the single particle basis. Put (f/Ar)nn' = A„5„„/. 
Then the above equation is verified by the following simple computation: 

N N 

n=l n=l 

N 

= nTr^i(-l)^^^^exp(/lAlnAj 

n=l 
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N 



= Tr^.(-irexp^/t/^lnA„. 

n=l 

The multivaluedness of the logarithm does not cause any ambiguity here as 
exp27ri/^/„ = 1. 

A similar expression for the inverse of a determinant can be obtained 
by substituting bosons for fermions. Let Bn be the bosonic Fock space 
generated by canonical boson operators 6]^, 6„ with commutation relations 
^n'^li ~ ^li^n' ~ ^nn'- For I7I < 1 wc havc the identity 

N 

Det(l-7C/;v)"' = TVB^exp bi{]njUj,)n„'bn'. 

n,n'=l 

Again, this is verified by diagonalizing Un and computing a single boson 

trace, 

Trs^e^i^'i'''^ - e°^^"^ + e^^^°^ + • • • + e"'"^ + • • • = (1 - A)"^ . 

The inequality I7I < 1 is needed in order for the sum on the right hand side 
to converge. 

Let us now assemble the various pieces to build a formula for fl. There 
are two determinants in both the numerator and the denominator, so we 
introduce two fermionic and two bosonic Fock spaces and B^, and write 

Det(l-7+i^Ar)Det(l-7-i^j^^) 
Det(l-7+of/iv)Dct(l-7_o?7^i) 

= Tre+^^+^g-^^-(-l)^^ exp^(6^„(ln7+o^jv)nn'^+n' 

nn' 

nn'b—n' 

+f+nO-^^+lUN)nn'f+n' 
+fU^n^-lU^\n'f-n') ■ (2) 

Here Np = J2n=i if+nf+n + f-nf-n] ; ^i^d Convergence of the sums requires 



|7-to| < 1- To obtain Q wc need to replace Un hj e^^U^ and integrate over 
6. If A"-!- = (^±n^±n + /±n/±n) couuts the uumbcr of retarded (+) or 
advanced (-) particles, this produces an extra factor 

^7 rift 
J 27r 
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under the trace. The operator P gives unity when apphed to a state I'lp) 
with (A^+ — A^_) \ip) = 0, and zero otherwise, and thus projects the super 
Fock space Bj^ J-'^ ® B]^ ® onto the subspace with equal retarded and 
advanced particle numbers. Hence, on writing J2nn' b+n{^^l+o^N)nn'^+n' —'■ 
h\\n (7_,_ot/7v)^+ simplify the notation, we obtain 

^(7+0, 7+1; 7-0, 7-i) = 

Tr^+=iv_(-l)'^"exp(6tln(7+ot^^)&4- + hl\n{^_,U^^)h_ 

+/iln(7+if/^)/+ + /lln(7_if/^^)/-) . (3) 

The final step in the construction of the supersymmetric formula for O 
is to pass from the trace over super Fock space to a Berezin integral, by the 
introduction of generalized coherent states. This step is perhaps unfamiliar 
and deserves a certain amount of explanation. 

Let me first illustrate the basic idea at the example of 7+0 = 7-0 = 0, in 
which case the bosonic degrees of freedom are absent. Thus, we start from 
the relation 

fi(0, 7+1; 0, 7_i) = Tr^+«^-^F exp [A In (7+if/^)/+ + /I In {l^,U^')f-) 

where Pp = {^-k)-^ ^^'^ dO ey^^ieY^nU+nf+n - f-nf-n)- Consider then the 
BCS (or pair) coherent states 

|Z)=exp( /!n^n„' /In' I |0) 

\n,n'=l I 

with complex amplitudes Z^n'- Being made from pairs of one retarded and 
one advanced particle, these states satisfy the constraint Pf\Z) = \Z). More- 
over, they enjoy the key property of providing a resolution of the projector 
Pf: 

Pf= J dfi^{Z,Z)\Z){Z\ , 

where dfXNiZ,Z) = const x Det(l + ZtZ)-^^-! n^„'=i d^nn' dZ„„/. The 
proof proceeds via Schur's lemma, by demonstrating that / dn^^Z, Z) \Z) {Z\ 
commutes with all of the operators f^nf-n' f^nf+n' {n,n' = 1, . . . ,N). 
A standard computation on coherent states yields 

{Z\ exp Y {ft In (7+if/^)/+ + /I In {l-iU^')f-) \Z) 

nn' 

= Det(l + 7+i7-i^^t/iv^t/iv')- 
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Hence, by inserting the coherent state resolution of Pp into the expression 
for Q{0, 7+1; 0, 7_i), we arrive at the formula 

1^(0, 7+1; 0, 7-1) = / dfij,{Z,Z)Bet{l + ^+,^_,Z^Uj,ZU^') . 

£NxN 

Our aim now is to extend this formula to allow for a nonvanishing first 
and third argument of Q, in which case the Fock space expression involves 
both fermions and bosons. This requires a supergeneralization of the notion 
of coherent states and their integration measure. The purpose of the next 
section is to review some of the basic mathematics needed. We will return 
to the task of expressing f2 as a coherent state integral in Sect. 

5 Basic notions of analysis on supermanifolds 

Familiarity with the elementary theory of differentiable manifolds is assumed. 
In the present section the summation convention is in force. 

Let AE"^ be the Grassmann algebra over M with q generators ip'^yip'^, . . . , ip'^. 
By definition, the Grassmann generators anticommute: ip^ip^ + ip^ip'^ = 0. 
They are known to physicists from their use in the construction of path 
integrals for fermions. Note in particular = —ip^ip^ = 0. 

Now pick some domain U G MP. A map F from U into the Grassmann 
algebra, 

F : U ^ AM."^ 

a ^ Fo(a) + F,(a)# + F,,,,(a)#lV'^^ + ■■■ + F,,...,,(a)^*^■■■^*^ 

will be referred to as a superfunction (or simply a function) on U. Fq is called 
the number part of F. We distinguish between even and odd superfunctions 
(the former are even and the latter are odd in the Grassmann generators). 
This distinction endows the algebra of superfunctions with a Z2-grading. 

We can now describe what is meant by a {p, g) -dimensional real-analytic 
supermanifold in the sense of Berezin, Kostant and Leites @| . Take a (real- 
analytic) manifold M of dimension p and cover it by domains Ui, U2, ■ ■ ■■ Let 
p even and q odd coordinate superfunctions be given on each domain. [For 
the purpose of illustration, we denote these by 

x\...,xP; e,---,e onf/i, 
and by . . . , y^; ?7\ . . . , r/'? on f/2 .] 
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Assume that transition functions exist on overlapping domains and are con- 
sistent with the Z2-grading. [For example, on Ui f] U2 this means that there 
are relations 

(and the corresponding inverse relations) with /q, /^;, ip\, film ^tc. being 
functions on a subset of W.] If the transition functions are analytic functions, 
we call the algebra of superfunctions on M generated by the coordinate 
superfunctions a real-analytic supermanifold of dimension {p,q). Complex- 
analytic supermanifolds are defined in an analogous manner (just replace 
Am'' by AC, the Grassmann algebra over C with q generators, and require 
the transition functions to be holomorphic). The manifold M is called the 
base of the supermanifold. 

Of course, a proper definition avoids any reference to specific coordinates, 
and in a more mathematical exposition one would define a supermanifold 
to be a "sheaf of algebras of superfunctions However, the coordinate 
description given here is good enough for our purposes. 

The standard example of a complex-analytic supermanifold is CP^'^. Let 
z^, be canonical coordinates of C^, and consider the algebra of superfunc- 
tions: C^\{0} AC^, with the complex Grassmann generator of Ac^ denoted 
by C} . Then focus on the subalgebra of superfunctions which are homoge- 
neous of degree zero in these generators, i.e. are invariant under rescalings 

—>■ cz^, — > C2^, cC} with c G C\{0}. This algebra is what is called 

the complex projective superspace CP^'^. Its base is the ordinary complex 
projective space CP^. Because the null element has been removed from C^, 
at least one of the two coordinates z^, z^ is nonzero. Let \J\ be the domain 
where z^ 7^ 0. Then homogeneity means 

f{z\z\C) = f{\,z'lz\C}lz'), 

so that / is a function of two variables z := z^ j z^ and C := C^/-^^^) which are 
taken to be the coordinate (super)functions on \J\. Similarly on [/2, defined 
by 7^ 0, we have 

/(^^^^c') = /(^v^^l,cv^'), 

and here we put z' := z^ j and C' '■= C^/z^- The transition functions on 
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C/i n are then given by 



z z 

The domains \J\ and Ui cover the base CP^, and the transition functions 
are holomorphic, so CP^'^ is a complex-analytic supermanifold, the complex 
dimension being (1,1). 

Of the many structures that exist on supermanifolds, the most important 
one for our purposes is the Berezin integral. To define it, we must first 
introduce the notion of Berezin form, also called superintegration measure 
or integral superform. If A denotes the algebra of superfunctions of our 
supermanifold and Q*M the space of volume forms on the base M, a Berezin 
form a; is a linear and local map 

u: A ^ n*M, 
f ^ uj[f], 

i.e. a linear and local rule for converting superfunctions into volume forms 
on M. This process of conversion is sometimes called the "Fermi integral" 
or "integration over the Grassmann variables" . 

For example, an interesting class of Berezin forms on CP^'^ is 

= d^ A dz o _ o (1 + + CO 

for n E N. The meaning hereof is this: to compute c<j'^"')[/], one first multiplies 
the superfunction / by {l+zz+(()^^\ then one takes two derivatives d'^/d(d( 
(thereby removing the Grassmann generators), and finally one multiplies by 
dz A dz to obtain a two-form on CP^. To illustrate, we take 

/ = /o(^,^) + /oi(^,^)C + /io(^,^)C + /ii(^,^)CC, 

in which case we get 

^(n)[^] = ( MhiL _ n /°°^^"M dz A dz . 

A special role is played by the Berezin form uj^^\ Converting cu^^^ from one 
chart to another by using the transition functions z — 1/z' and ( — C'/^', one 
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finds that cu^^^ operates by the same expression in both charts. (Moreover, 
cu(i) is invariant under an action of the Lie supergroup SU(2|1) on CP^'^.) 

The definition of Berezin forms immediately leads to a notion of integra- 
tion on supermanifolds called the Berezin integral, or superintegral. Since a 
Berezin form u converts a superfunction / into a volume form a;[/] on M, the 
natural thing to do is to integrate a := uj[f] in the ordinary fashion to pro- 
duce the number a. Thus a Berezin integral is defined to be a composite 
map 



JM 

taking sup erf unctions into the complex numbers. Note that this definition 
decrees the Berezin integral to be a two-step process: first the "integra- 
tion" (or, rather, differentiation) of the Grassmann generators is done, and 
afterwards the differential form uj[f] is integrated in the usual sense. The 
advanced user occasionally finds it convenient to deviate from this rule and 
perform part of the ordinary integrals first. However, in case of doubt the 
definition recited above is the one to go back to. 

Ordinary integration is coordinate independent, i.e. one has the option 
of changing integration variables by the substitution rule. It turns out that 
this option also exists for superintegrals, albeit with a certain caveat. To 
describe the relevant statement, which is due to Berezin, we introduce the 
short hand notation 

D{x,0 :=dx^A...Adx-A---A, 

and consider a change of variables (consistent with the Z2-grading) from x, ^ 
to y, T] by 

= ^^{x\...,xP-e,...,^') (j = l,...,g). 
The role of the Jacobian of ordinary analysis is taken by the Berezinian, 

-l^j :=SDet(^^^,^^^^. 

where the superdeterminant SDet of a supermatrix ^ ^ ^ ^ (with the ma- 
trix elements of A, D being even and those of B,C odd) is defined by 
A B \ Bet{A - BD-^C) DetA 



SDet 



CD) DetD Det{D -CA-^B) 
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Berezin's theorem fOl states that the substitution rule is vahd, i.e. 



provided that / is compactly supported. 

The condition of compact support, unfamiliar from ordinary analysis, 



invites some explanation. For an instructive example |]I0|, consider the non- 
compact supermanifold of functions: ]0, 1[— AR^, with canonical coordinates 
X, If we choose to integrate with the Berezin form dx d"^ / d^^d^^ , a su- 

perfunction / that depends only on the combination x + ^^^^ has the Berezin 
integral 



1 



On the other hand, if we made a change of integration variables 
which is easily seen to have unit Berezinian, we would get 







Because the Grassmann derivatives d"^ /drfdr]^ have nothing to act on, the 
right hand side vanishes although for equality with the left hand side, it 
ought to be /(I) — /(O). However, both sides vanish, and the discrep- 
ancy disappears, if we take / to be compactly supported (which implies 
/(O) = /(I) = 0), as is required in order for Berezin's theorem to apply. An 
alternative way of fixing the problem is to pass from the noncompact interval 
]0, 1[ to the compact circle S^, by imposing periodic boundary conditions on 
X, so that /(I) - /(O) = 0. 

The above example signals a general complication, which one has to con- 
front when changing variables in a superintegral. While ordinary volume 
forms transform simply by the Jacobian, the transformation law for Berezin 
forms is not given only by the Berezinian, but involves an additional, "anoma- 
lous" term /3: 

D{y,^)=D{x,i) Ber (^) + , 
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whenever some of the even coordinates are shifted by nilpotents. In our 
simple example, 



Diy,v)f 



D{x,Of + dx^f 



The anomaly f3 always has the property of being exact (in the sense of dif- 
ferential forms), i.e. p[f] = da[f] where da[/] is the exterior derivative of 
some [p — l)-form a[f]. Therefore, by Stokes' theorem the anomaly can be 
integrated and converted into an integral over the boundary of the base of 
the supermanifold (or the boundary of the chart used, if there exist coor- 
dinate singularities). Rothstein |jlO[ has given an explicit expression for the 
anomaly in terms of the vector field generating the transformation from x, ^ 
to y,r]. 



6 Supersymmetric formula for maps (II) 

After this brief tour of some basic supermathematics, we return to the task 
of setting up a supersymmetric formula for the pair correlations of quantized 
symplectic maps. We left off at Eq. (^ expressing f2 as a trace over S, the 
subspace of super Fock space B'li®T^®Bji®J^'^ determined by the condition 

= (equal number of retarded and advanced particles). 

As was motivated at the example of the purely fermionic case in Sect. ^, 
our goal is now to trade the trace over S for an integral over generalized 
coherent states. To that end, we observe that every Fock state I?/') satisfying 
(A^+ — A^_) \ip) = can be obtained by repeatedly acting on the vacuum with 
the pair creation operators &U&1„,, /+n&L„', and flnfln' i^^^' = 

1, . . . , N). This leads us to consider coherent states of the pairing form 

\Z) := exp Y: {blnZT'bln' + blnZ^fL' + flnZl^'hl^, + flnZZfl^) |0) , 

nn' 

where mathematical consistency requires taking complex numbers for the 
matrix elements Z^^,, Z^^, and Grassmann generators for Z^^,, Z^^,. It is 
convenient to assemble these into a supermatrix 
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Then a short hand notation for the pair coherent states is 

\Z) =exp(4ZcM|0) 



where the creation operators c± stand for both foj. and /±. We will also need 
the dual of a pair coherent state 



exp (c_ctZc+), Z=( ^ff 1 ' ^ 



^FB ^FF j ' ^ - I -Ijv 



In |0] it was shown that, with a suitable choice of Berezin form D{Z^ Z), 



the pair coherent states provide a resolution of the projector P, 
P = j -^exp^^(iV+-iV_)= jD{Z,Z) {Z\Z)-^\Z){Z\ 



where the domain of integration is fixed by the conditions 

^FF^_^FFt^ ^BB^^^BBt^ 1 - Z^^t^BB > q _ 

This resolution follows from Schur's lemma, once it has been 

demonstrated that the right hand side commutes with all bilinears in Fock 
operators that commute with the projector P. The crucial relation needed 
for the latter statement is the invariance of the Berezin integral under a left 

translation through g = i ^ D ) ^ GL(2A^|2A^), 



D{Z,Z)f\Z,Z) = 
j D{Z, Z)f[{AZ + B){CZ + D)-\ {C + DZ){A + BZ)-^] 

which determines the Berezin form D{Z, Z) to be the invariant superintegra- 



tion measure on Efetov's model space II (or, using the terminology of |12|, 
on a Riemannian symmetric superspace of type AIII/AIII). The invariant 
Berezin form D{Z, Z) turns out to be flat as a result of cancellations due to 
supersymmetry: 

D{Z, Z) = Yl dZ^^,dZl^,dZll,dZll, + ... 
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modulo boundary anomalies, which are indicated by the dots. The normal- 
ization oiD{Z,Z) is fixed by 1 = (0|P|0) = J D{Z, Z){Z\Z)-\ 

Given the coherent state resolution of the projector P, the trace of any 
operator O on Fock space can be converted into a Berezin integral: 

Tr^+=jv_(-l)'^"0 = STrPO = J D{Z, Z){Z\0\Z) / {Z\Z) . 

Application to Eq. (^), followed by a computation of {Z\Z) and {Z\ . . .\Z), 
results in the desired formula for VL: 

(](7^;7_) = j D{Z,Z) SDet(l-ZZ)SDet-i(l-Z(7+®f/^)Z(7_®f/^i)) . 

(4) 

Here 7± = diag(7±05 7±i) are diagonal 2x2 matrices, and 7+ ® Un is meant 
to be a supermatrix: 



7+ ® Un 



7+of/Ar 

1+iUn 



A closely related derivation of the formula for is provided by the so- 
called color-flavor transformation (|jll| and The two-level correlation 
function R2 follows from Vt by taking two derivatives, as before. 



7 Semiclassical limit? 

So far, we have achieved no more (and no less) than an exact reformulation 
of the original problem. As should be clear from its derivation, the super- 
symmetric formula for the correlator VL applies to any unitary matrix U^, 
independently of whether this matrix arises from quantization of a map or 
not. In order to go further and extract nonempty information from our for- 
malism, we need to exploit the fact that Un has a semiclassical limit and, in 
particular, we will have to make a distinction between chaotic and integrable 
systems. 

To prepare these steps, we write VL = J D{Z, Z)e^^ where, borrowing the 
terminology from quantum fleld theory, the function 

S = -STr ln(l - ZZ) + STrln (l - ^(7+ ® f/jv)^(7- ® Un^)) (5) 

will be called the "action functional". Now recall that the BB- and FF- 
blocks of Z are complex N x N matrices. Viewing them as linear operators 
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on the Hilbert space 'Hn, we might say that they are similar to the den- 
sity matrix of a quantum mechanical state. And, in fact, an expression like 
Z^B' := UnZ^^U^^, which occurs under the second logarithm in S*, can be 
interpreted as being the "density matrix" evolved by one (inverse) time 
step. One may therefore be tempted |Q to transform to a Wigner represen- 
tation, postulate a semiclassical limit for Z and subject the action functional 
to a semiclassical expansion. However, such an expansion is entirely un- 
controlled in the present context and is, in fact, false. Unfortunately, a lot 
of confusion has been created among the community (including myself) con- 
cerning this point. Let me therefore change the style of presentation and give 
a detailed exposition of the issues as I see them. The reader is warned that, 
given the current level of understanding, the following has to be somewhat 
qualitative. 

If A is a linear operator (on some quantum mechanical Hilbert space) with 
kernel A[q, q') in, say, the position representation, one defines the Wigner 
transform (or Weyl symbol) of A by 

^A(g, P) = E Mq + 972, g - q'/2) e^^^'^ 

whenever the transform exists. (In our case there is a technical complication 
due to the fact that the "position" on a compact phase space is not glob- 
ally defined, but this is a peripheral issue and we are not going to worry 
about it here.) Now if A and B are pseudodifferential operators with Wigner 
transforms OAifl^V) and crB{q,p), their operator product has a semiclassical 
expansion 

(TAB = (TAfTB + ifl/2t){aA, CTb} + C(?l'), (6) 

where ctactb is a product of functions, and the curly brackets denote the 
Poisson bracket. 

The question now is whether this kind of expansion applies to products of 
the supermatrix Z. To be sure, there is nothing that prevents us from trans- 
forming Z to a Wigner-like representation, but are we really allowed to make 
a semiclassical expansion of the action functional?! The question is a valid 
one, for the semiclassical expansion (^ does require A and B to be pseu- 
dodifferential operators. Without the smoothness property provided by this 
condition, the Wigner transform of the operator product AB does not sepa- 
rate into the product of Wigner transforms a AO's in the limit h ^ 0. For the 
standard operators appearing in quantum mechanics (position, momentum. 
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energy etc.) the smoothness condition is of course always satisfied. However, 
there is no principle that guarantees smoothness of the Wigner transform az 
in the present case. Indeed, Z is not a fixed observable but is a variable of 
integration. While a small fraction of the Z that are integrated over do cor- 
respond to smooth az, the vast majority do not. (Recall that the integral is 
over all complex supermatrices Z, Z with the integration measure being that 
of a Riemanian symmetric superspace.) This invalidates the error estimate in 
d^). One may still cherish the hope that the dominant contributions to the 
Z-integral do come from smooth configurations. Admittedly, this is what 
happens for Feynman's path integral in the limit h ^ 0, which is dominated 
by paths that are extrema of the action functional. However, the present 
situation is quite different. There exists no saddle-point or stationary-phase 
or other principle here that would favor smoothness in general. 

To recognize the severity of the problem, recall that the unitary operator 
Un has eigenvalues e*^" and eigenfunctions ipn, and there exist N of them. 
As a result, the conjugation Ad{U]\i) : p t-^ Uj^pU^^ has eigendensities 
p = ipn ® "fpn with eigenvalues e**^^""^"-' = 1. To compute the pair corre- 
lation function in the microlocal scaling limit, we need to set 7± = e''^^^^ 
and let A^ go to infinity. In this limit the action functional acquires A^ zero 
modes given hy Z = Y,n=i C^'^'^ipn^n, where the C^"^ are 2x2 supermatrices 
and are otherwise arbitrary. Thus there exist A^ directions through the point 
Z = Z = along which the integrand is exactly neutral - a circumstance that 
surely causes very-large-amplitude fiuctuations about this point. Moreover, 
the zero mode densities ipn'^Pn are never smooth. Indeed, in the integrable 
limit they are sharply localized on invariant sets, while in the chaotic limit 
they are highly irregular functions that vary on the shortest scales permitted 
by the uncertainty principle. In other words, the action functional @ is un- 
stable with respect to a large number A^ of zero modes and none of these has 
a smooth Wigner transform. In view of all this, we had better abandon the 
hope that the dominant contributions to the Z-integral come from configu- 
rations that possess the smoothness of a pseudodifferential operator. Quite 
on the contrary, the most important configurations are those where az varies 
with the shortest wave length possible! As h is lowered, the dominance of the 
A^ zero modes keeps introducing an increasing number of relevant field con- 
figurations that vary on finer and finer scales. In order for the semiclassical 
asymptotics (|]) to set in, we would have to intervene and limit the variations 
of az to wave lengths that are greater than some fixed (i.e. /i-independent) 
minimal scale, and then let h go to zero. Without such an intervention, the 
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Poisson bracket (and higher) terms in fail to be /i-independent but diverge 
(with increasing powers) as h —* 0, so that truncation of the semi classical ex- 
pansion of S is totally unjustified. Consequently, any conclusion resulting 
from the use of this expansion in (|]) is at risk to be false. 

Given the invalidity of the semiclassical step, how can we use the super- 
symmetric formula (^) to develop the theory further? Not surprisingly, the 
answer depends on what we are trying to achieve. Here we are pursuing no 
more than the modest goal (modest from a physicist's perspective) of formu- 
lating a precise version of the BGS conjecture and substantiating it. As will 
be shown in the sequel, this restricted goal offers us the option of postponing 
the semiclassical step until the very end of the calculation. 

8 Regularization 

Notwithstanding the fact that the formula @ is well-defined and, in fact, 
mathematically rigorous, we are facing a key difficulty: the existence of 
one-parameter groups along which the integrand lacks stability prevents us 
from investing semiclassical input and makes it difficult if not impossible to 
compute something from (^) as it stands. To make any progress at all, we 
must first solve this stability problem. Consider therefore, in the microlocal 
scaling limit 7± — > 1, the quadratic part of the action functional at Z = Z = 
0: 

^(2) = STr ZZ - STr ZUjvZU^^ = STr Z{1 - Ad(f/jv))Z , 

which identifies the Hessian at this point as the operator 1 — Ad(f/Ar). Now 
recall that the classical limit of Ad(f/Ar) is x*? the Frobenius-Perron operator 
of the inverse map To prepare the discussion of the quantum mechanical 
operator Ad{UN) in the limit A^ —>■ oo, we shall summarize what we know 
about X*- 

First of all, because x* is unitary with respect to the Liouville L^-measure 
on phase space, the spectrum of x* li^s on the unit circle in C. (If the map 
is mixing, the spectrum is known to be absolutely continuous |0.) Second, a 
salient feature of this spectrum is its instability with respect to regularization 
by smoothing. More precisely, if we project x* from the full L^-space to a 
subspace of smooth functions, by coarse graining the phase space or imposing 
an ultraviolet (UV) cutoff, unitarity is lost and the spectrum moves to the 
interior of the unit circle. Third, if the map x is mixing, the spectrum of the 
projected operator consists of one isolated nondegenerate eigenvalue at -|-1 
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(with the corresponding eigenfunction being the invariant density po = const) 



and so-called Ruelle resonances [|TJ] inside the unit circle. It is important 
that this feature persists, i.e. the unitary spectrum is not recovered, when the 
UV cutoff is lowered to zero. (This is a statement of non-interchangability 
of limits.) Note also the difference between the quantum and classical cases: 
the quantum operator 1 — Ad(f/7v) always has zero modes, whereas its 
classical limit 1 — x* on smooth functions only has a single zero mode for a 
mixing system, namely the invariant density. 

The usual approach taken in the study of Ruelle resonances is to regularize 
in the ultraviolet by imposing a smoothness or differentiability condition on 
the phase space functions acted upon by the Frobenius-Perron operator. A 
similar effect is achieved if we leave the functions unchanged and, instead, 
implement the UV regularization on the Frobenius-Perron operator itself. 
For our purposes the latter procedure is more convenient. What we do is to 
make the replacement 

X — ' r X , 

where the operator r*^ is chosen to be the exponential of a sum of squares of 
Hamiltonian vector fields S^, 

s 

r" = exp e ^ . 

k=l 

For example, for M = we take for the generators of rotations {k = 
1, 2, 3), in which case S| is simply the Laplacian on S^. More generally, 
the choice of Hamiltonian vector fields is constrained by the requirement that 
Y^k'^k t>e an elliptic operator. The effect of r^, then, is to regularize in the 
ultraviolet by suppressing fluctuations on short scales. It is important that 
r*^ is the exponential of a differential operator of second order, whereas the 
Frobenius-Perron operator can be thought of as the exponential of a first- 
order differential operator. Therefore, is always relevant in the ultraviolet, 
no matter how small is the value of e. As a result, remains effective as 
a regulator even in the limit e — > and, again, the unitary spectrum is not 
retrieved when the UV cutoff is removed. 

In the following we will make use of the fact that, if the map x is mixing, 
the operator 1 — r^x* is strictly positive on the orthogonal complement of 
the uniform state po- The physical reason for such behavior is this. Consider 
the long time dynamics obtained by iterating r'^x* ^^id, for simplicity of 
the argument, assume x to be uniformly hyperbolic, which means that x is 
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everywhere stretching along unstable manifolds and contracting along stable 
ones. Thanks to the eUipticity of J2k=i'^h the operator r*" counteracts the 
contraction along the stable manifolds via diffusion. After many iterations, 
the combination of stretching by x and smoothing by r'^ will attract any initial 
state to the uniform state. Thus, as time / goes to infinity, the operator 
^j,e^*y approaches the projector onto the uniform state po- This implies 
strict positivity of the real part of the spectrum of 1 — r^x* operating on the 
orthogonal complement of po- 

Let us now turn from the space of functions on M, acted upon by x*; 
to the quantum mechanical space TIn ® 'H*^^, acted upon by the operator 
Ad{U]\f). Crudely speaking, the quantum space differs from the classical one 
by a finite resolution, or cutoff, set by Planck's constant. When this cutoff is 
much smaller than the regularization scale set by e, we have a good semiclas- 
sical expansion, (^), and the regularized Ad(f/Ar) spectrum looks qualitatively 
similar to the regularized Frobenius-Perron spectrum. On the other hand, 
if h is kept fixed at a finite value, the destruction of the unitary spectrum 
of Ad(Uj\f) with increasing e is not abrupt but is rounded off. If we take 
e ~ with a > 2d, the strength of the perturbation caused by the regulator 
becomes smaller than the level spacing of Ad(f/Ar) [which is 2-11 /N"^ ~ h'^'^], 
and is therefore negligible, in the limit h ^ 0. In contrast, if we choose the 
scaling exponent a in the range < a < 1, then by power counting the 
regulator remains relevant for h ^ 0. This is good news: it suggests that 
we can kill the unitary spectrum of Ad(?7iv), and thereby get rid of the zero 
mode problem if the classical map is mixing, by a regulator with strength e 
that vanishes in the classical limit. 

The basic strategy to follow should now be clear: our aim must be to sub- 
stitute for the unitary operator Ad{UN) in the quadratic functional S^"^^ the 
regularized operator i?^Ad(?77v), where i?^ is a quantization of the diffusion 
operator r^ Such a substitution is expected to turn the unitary spectrum of 
Ad(f/Ar) into something akin to the Frobenius-Perron spectrum of the classi- 
cal operator x*- If X is mixing and e depends on h as e ^ with < a < 1, 
only the invariant zero mode of 1 — Ad(f/Ar) remains neutral. The remaining 

— 1 zero modes turn into positive directions of the regularized Hessian 
and therefore become amenable to perturbative treatment. From experience 
with the application of the supersymmetric method to disordered electron 
systems we then anticipate that a careful integration over the invariant 
zero mode will give the random matrix answer for the correlation functions. 

What is unsatisfactory about this strategy is the absence of any explana- 
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tion of where the UV regularization came from. So far, regularization seems 
to be an ad hoc procedure, or "dirty trick", introduced so as to yield the 
result one wants to get. The situation would be different if we were working 
in a quantum field theoretic context. In that case, UV regularization would 
be well justified by the fundamental and inevitable existence of ultraviolet 
singularities (signaling incompleteness of the field theory, or "new physics" 
at short scales), which must be tamed by the introduction of a cutoff. In 
contrast, our formula has no singularities whatsoever but is perfectly 
well-defined - we just don't know how to evaluate it - and there is no a 
priori need for regularization. In passing we mention that, if we ignored the 
arguments of Sect. |^ forbidding us to apply the semi classical expansion (||) to 
the action functional S, and applied the expansion anyway, then we would, 
in fact, end up with an ultraviolet singular field theory [Q. However, the UV 
divergencies so introduced are entirely artificial - they simply confirm the 
fact that we have made improper use of the semiclassical expansion - and do 
not justify the postulate of any regulator. 

To reiterate, it is illegal just to postulate UV regularization as a technical 
device or trick. For a convincing argument, one needs to explain exactly 
what modification of the original problem is implied by regularization, and 
one has to do that in a frame that precedes the supersymmetric quantum 
field theory formalism. As will be argued in the remainder of this section, 
regularization of the Hessian of (^) at Z = Z = amounts to averaging over 
the "quantum" ensemble described in Sect. ^. 

Recall that, if {k = 1, s) are uncorrelated Gaussian random variables 
with zero mean and variance (^|) = 2e, the ensemble is specified by 



By averaging the Hessian 1 — Ad(f/jv(0) over the probability space of the 
parameters ^k, we get 






exp iJ2^kad{Xk,N)/fi) =exp -e^ ad^(Xfc,7v)/ft^ + ••■ ■ 



k I k 
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[Note that, since the integrand in depends on 11^ in a nonhnear way, 
later we will have to confront averages of powers of Ad {U]sf{^)). For the time 
being, we concentrate on the Hessian.] The terms indicated by dots in the 
last expression become negligible in the semiclassical limit h ^ 0. Taking the 
commutator with Xk^N in this limit corresponds to applying the Hamiltonian 
vector field H^: 

ad(Xfc_7v) = [Xk,N, •] — ^ ifi^k , 

r% as desired. 

To conclude, averaging over the specified ensemble regularizes the Hessian 
in the ultraviolet. As was argued earlier, if the classical map x is mixing the 
real part of 1 — r'^x* is positive on all states that are orthogonal to the uniform 
state pq. In view of the semiclassical limits i?^ — > and Ad(f/Ar) — > x*, we 
expect the same to be true for the regularized Hessian 1 — i?^Ad (Un) on 
Vj^, the subspace of TCn^'H*^^ orthogonal to the invariant element. This key 
observation allows us to proceed as follows. 



9 Asymptotic expansion 

We set 7± = e^""^*/^ where c± = diag(c±05 c±i) are 2x2 matrices and, to 
abbreviate the notation, we put Zu := {l2®Uj\f)Z {12^17^-^). Then, since c± 
are held fixed while N is sent to infinity, we may expand the action functional: 



S = -STr ln(l - ZZ) + STr ln(l - ZZi 



u) 

--^STr {c+ZuZ{l - ZuZ)-' - c.ZZu{l - ZZuY^) + 0{llN^) . 
To isolate the invariant zero mode, we make the substitution 

Z = (C + Zo)(l + ^oC)"\ 

z = (c + lo)(i + ^oC)-'- 

The supermatrices C, and C, are taken to be traceless in each block (i.e. 
Tr = Tr C^"'^ = for a, r = B, F), and the zero mode is parameterized by 
Zq = zq ^ In, Zq = zq ^ In with 2x2 matrices 

/ ^BB ^BF \ / ~BB ~BF 

_ Zq \ ~ _ / ^0 ^0 

^0 — I ^FB ^FF ) 5 ^0—1 ~FB ~FF 
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Because the transformation 

is an isometry of the Riemannian symmetric superspace, the Berezinian of 
the change of variables from Z to (C, Zq) is unity and -D(Z, Z) factors into a 
product of Berezin forms D{zq, zo)D{(, Q. (We pay no attention to possible 
boundary anomalies here.) Now, by inserting the expansion of the action S 
into (^), changing variables as indicated, and using identities such as 

(1 - ZuZ)-' = (1 + ZoC)(l - (uO-\l + (uZo)il - ZoZo)-' , 
we obtain 

^^^^c+/N. ^-^c^/N^ ^ J ^(^^^ n,ic+/N, C./N; Zo, Zq) 

X exp iSTrcm (c+Zo5o(l - zoZq)'^ - c^zqZq{1 - ZqZqY^^ 

with Oi defined by 

^]l(c+/iV, c^/N-z^rzo) 
= J Die, SDet(l - CC)SDet-i(l - CCf/) 

X exp (^STr c+(l + ZoC)(l - Cf/C)"'Cf/(C + ^o)(l - ^o^o)"' 

- ^STr c_(l + ZoCu)il - CCu)-'aCu + ^)(1 - ZoZ^)-') . 

Note that convergence of the last integral is ensured by the condition Imc+o > 
> Imc_o resulting from |7±o| < 1- 

Our next goal is to establish an asymptotic expansion for the function 
Qi. To that end, we temporarily set c+/N = c_/N = and consider the 
integral 

ni(0,0;zo,5o) = y'z^(C,C)SDet(l-CC)SDet-i(l-CCc/) • 

To compute this, we observe that both the Berezin form D{(,() ^i-nd the two 
superdeterminants are invariant under transformations 
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with A = a ® In and D = d ® 1^. (Indeed, such transformations com- 
mute with Ad(f/Ar) and are isometrics of the symmetric superspace.) This 
invariance resuhs in the value of the integral being equal to the value of the 
integrand at the origin ( = ( = 0: 

Qi{0,0; zo,zo) = 1 , 

by what is often called the Parisi-Sourlas-Efetov-Wegner theorem in dis- 
ordered electron physics. The mechanism underlying the theorem is called 
"localization" in the mathematics literature. In a setting closely related to 
equivariant cohomology, localization of superintegrals has recently been dis- 
cussed by Schwartz and Zaboronsky |T^. With reference to their results we 
can argue as follows. On setting 

A = expl Q 1 (g) Itv , D = exp ( ^ ^j^l^, 

and differentiating with respect to the Grassmann parameters a, a, 5, 5, the 
transformation ( i— ^ A(D~^, ( i— >■ D(A~^ gives rise to four odd vector fields, 
each of which leaves the integrand invariant. The joint zero locus of these 
vector fields is ( = ( = 0, and the superdeterminant of the Hessian at this 
point is unity by cancellation due to supersymmetry. Moreover, the vector 
fields are compact in the sense of Schwartz and Zaboronsky. Therefore, by 
the version of the Parisi-Sourlas-Efetov-Wegner theorem proved in |]T5|, the 
integral equals unity as claimed. 

Given the result ^1(0, 0; zq, zq) = I, one would like to proceed and make 
a saddle point approximation around ( = ( = 0, to produce an asymptotic 
expansion for Qi{c+/N,c_/N] Zq, Zq) of the form 1 + A^~^/(c+, c_; Zq? ^0) + 
0{N~^). However, this is impossible as it stands. The reason is the notorious 
problem of — 1 zero modes: there exists an {N — l)-dimensional maxi- 
mal torus of group actions that commute with Ad(f/Ar) and leave SDet(l — 
CC)SDet^^(l — CCu) invariant, thereby causing a lack of stability of the can- 
didate saddle point ( = ( = 0. Fortunately, given all the preparations that 
were made in Sect. ^ we immediately know how to fix this stability problem. 

We substitute Un{C,) for Un and average over the probability space of 
the ensemble, replacing Q by (Q), and fli by (fli). The integral for Qi 
converges and, as a matter of fact, converges uniformly in the parameters 
C,k, if Imc+o > > Imc_o- Therefore, we may interchange the order of 
integration and ensemble averaging. By the same localization argument as 



26 



before, 

(fii(0, 0; ^0, z,)) = I Die, (SDet(l - CC)SDet-i(l - CCu^^))) = 1 • 

What has improved in comparison with the situation without ensemble av- 
eraging is that now we do have a good asymptotic expansion. The reason is 
that integrals such as 

I Dice) (SDet(l -CC)SDet-i(l -CCl/^.(o) 
xSTr c+(l - Ci/.(oC)-'Cc/.(oC(l - ZoZo)-') 

are expected to exist and remain finite in the limit ^ oo. This can be 
checked by expanding around the saddle point = = in the usual manner. 
By doing the resulting Gaussian integrals over (, (, one obtains traces involv- 
ing powers of the inverse of the regularized Hessian 1 — Rj^AdiU^) acting 
on Vj^. We have argued earlier that the regularized Hessian remains strictly 
positive in the semiclassical limit if the classical map is mixing. Therefore, 
the inverse of 1 — -R^Ad(f/Ar) exists on Vj^ and is uniformly bounded in A^. 
Moreover, as follows from semiclassical trace formulas, the traces encoun- 
tered are bounded uniformly in if the fixed points of the classical map are 
isolated. 

What all this amounts to is that, unlike fli, the ensemble average {Qi) 
does have an asymptotic expansion: 

mc+/N, c./N; zo, ~zo)) = 1 + N~'fic+, c_; Zo, So) + OiN~^) . 

Using this in the formula for the correlator (Q) we obtain 

lim (fi(e''=+/^;e-*"-/^)\ 

= / i„) exp iST. (e,.„i„(l - ..S.)- - c.S„.„(l - i„.„)-) . 

Computation of this definite superintegral is a standard exercise and gives 
the random matrix answer (|l]). The result for the pair correlation function 
follows by taking two derivatives. The dependence on Zo,Zq of the 1/A^ 
correction to can be shown to be such that the integral over zq, zq exists 
and the coefficient of 1/A^ in {fl) is finite. 

In conclusion, I believe that the arguments presented here are strong 
evidence in favor of the precise version of the BGS conjecture proposed. To 
turn them into a proof that is rigorous by mathematical standards, three 
improvements are required. 
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1. One needs to do a spectral analysis of the operator Rj^Ad^Uj^f) to sup- 
port our intuition of what happens in the limit — oo. 

2. The transformation of Berezin forms D{Z, Z) D{zo, zo)D{(, () needs 
to be scrutinized carefully, to rule out the possible existence of bound- 
ary anomalies (Sect. ^ that might interfere with our argument. 

3. A complete computation of the corrections is called for. In par- 
ticular, one must resolve the puzzle why the calculation of Andreev 
et al, which is closely related to ours, fails to reproduce the diagonal 
approximation to the double sum over periodic orbits. 
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Systems and Quantum Chaos" at the Newton Institute, Cambridge, where 
these lectures were delivered, I enjoyed a number of illuminating discussions 
with Steve Zelditch, who got me straightened out on the subject of semiclas- 
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Hamiltonian flow that averages to a diffusion operator. 
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